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ABSTRACT 


The  interface  between  a  solid  and  a  liquid  is  modelled  by  a  flat 
surface  with  an  array  of  sticky  sites  which  could  be  placed  on  a  regular 
lattice,  or  also  randomly.lt  is  first  shown  that  the  thermodynamics  and 
distribution  functions  can  be  expressed  entirely  in  terms  of  the 
distribution  functions  of  the  system  without  the  sticky 
sites . Furthermore ,  the  problem  of  the  occupation  of  the  sites 
(equivalent  to  the  adsorption  isotherm) , is  a  lattice  problem, which 
exhibits  phase  transitions  and  critical  points. A  simple  application  to 
the  interface  between  a  solid  and  a  fluid  of  hard  spheres  is  discussed. 


INTRODUCTION 


Real  fluid-solid  interfaces  have  a  good  deal  of  structure  and 
therefore  are  very  difficult  to  treat  in  Statistical  Mechanics. 
Traditionally,  the  literature  has  dealt  with  the  stucture  of  the  solid 
side  ,  and  therefore , ignored  the  liquid, or  else, treated  the  liquid  using 
accurate  integral  equations .but  assuming  that  the  surface  is  perfectly 
smooth. In  the  presept  work  we  discuss  a  model  of  the  interface  in  which 
the  fluid  is  treated  using  an  integral  equation, but  instead  of  a  smooth 
interface, we  have  an  array  of  sticky  sites  which  represent  the  atoms  in 
the  solid  side  wf  the  interface . This  model  is  still  tractable,  and 
incorporates  the  surface  structure  and  its  effects  into  the  theory. 

In  previous  work  (1]  an  exact  solution  of  this  model  for  the 
2  -  dimens ional  one  component  plasma  was  presented. 


The  model  of  Langmuir  [2]  was  perhaps  the  first  model  of  this  kind 


^  . 

V.\ 


buc  in  which  the  fluid  was  treated  as  an  ideal  gas. Further  treatments 

incorporated  the  structure  of  the  surface  as  a  perturbation  [3, 4], The 
adsorption  on  a  uniformly  sticky  surface  was  discussed  in  the  pioneering 
work  of  Perrara  and  Smith  [5]. In  our  model, the  sticky  potential  of  Baxter  u^(r) 
[6]  is  used. Consider  for  simplicity  a  fluid  of  spheres  witl^a  hard  core 
of  diameter  tf  and  density  j>  near  the  wall. The  position  of  the  center 
of  the  sphere  r-(x,y, z) .where  z  is  the  distance  to  the  surface,  and 
R-(x,y)  is  the  position  of  the  projection  of  the  surface. The  sticky 
interaction  is  then  given  by 


I 

fc 


exp[ -^u  (r)  ]-  1+x  (R)  o(z-<f 


/2) 


(1) 


where  X  (R)  is  the  stickiness  parameter .which  could  be  an  arbitrary 
function  of  R.In  most  cases  it  will  be  periodic,  and  in  the  case  in 
which  the  adsorption  sites  are  points  we  have 


A(R)  -  X  J  S (R  -n  a  -n  a 


(2) 


where  n^ .n^are  natural  numbers  and  a^.a^are  the  lattice  vectors  of  the 

sites  on  the  surface. 

* 

Consider  now  a  fluid  of  N  particles  in  a  container  of  volume  V. 
The  Hamiltonian  of  the  system  is 


H  -  H  +  H, 


(3) 


where  H  is  the  Hamiltonian  in  the  absence  of  the  sticky  interactions 


and  che  sticky  interactions  are 


V/  V*> 


where  che  sum-  is  over  all  sices  of  che  surface. The  canonical  parcicion 


function  is 


Zj:  (1/N!)  j  exp  (-0H0)  n  {(1+  A(R£  6  (z  . -<S/2) )  dr.] 

j  i-s., 


and  afcer  some  simple  manipulacions  we  get  [1] 


Z  -2°  J  (  //s!)p*(R..  R)\ 

rf  /  5  <-»  • 
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where  we  have  used  che  definition  of  the  n-particule  distribution 


function 


..s)  - 


)  "'/|TN's)  IZ.°Jj^dr^  .  .dr  expf-^H*) 


s)-  n(»#<i))  g  (l...  S) 

i~\  1  ^ 


It  is  clear  from  (5)  that  X  can  be  identified,  up  to  a  multiplicative 

constant  with  the  fugacity  of  the  adsorbed  particles . The  average  number  of  adsorption 
sites  <N  >  is  given  by 

<N  >-  In  Z  In  X  (9) 

3l 

Now, the  fraction  of  occupied  sites  is  ©  -N^/M  where  M  is  the 
number  of  sites  in  the  area  A  of  the  interface . The  number  of  occupied 
sites  per  unit  area  is 

n  -  <N  >/A  -  ©  w  (10) 

«.  a 

were  w-M/A  is  the  number  of  sites  per  unit  area. 

The  important  observation  is  that  the  partition  function  (6) 
is  isomorphic  to  the  grand  canonical  partition  function  of  a  lattice 
with  sites  distributed  on  the  surface  in  the  same  positions  as  the 
sticky  sites  of  our  model. The  interactions  amongst  the  particles  on  the 
sites  are  obtained  by  introducing  the  potentials  of  mean  force  of  the 
smooth  wall  problem  (SWP)Lf73 

gjd-..  s)-  exp  ( -0^(1,  j)  +  ^Tw(ijk)  +  w(ijkl)  +  .  .  ] ) 

'  (ii) 

where  the  potentials  of  mean  force  w(ij ) ,w(ijk) , . . . are  defined  in  terms 
of  the  SUP  interactions . The  indices  i,j,k,  refer  now  to  the  sites  on  the 
surface .  Furthermore  the  '.external,  field’  in  the  mapping  is  given  by 

a,  &  A*.  -  -  8-^  l  > 

(12) 

where  p°  (R)  is  the  density  at  a  given  site  R  on  the  surface. 


We  will  be  concerned  with  the  occupation  number  of  the  sites  on  the 
surface .There  are  many  different  systems  of  physical  interest  that 
can  be  studied,  but  for  the  sake  of  simplicity,  we  will  first  assume 
that  the  adsorbate  molecules  of  the  fluid  are  smaller  than  the  site 
separation. That  means  that  one  molecule  cannot  occupy  two  sites  at  the 
same  time,  and  therefore  our  problem  can  be  mapped  into  an  Ising  problem 
with  arbitrary  interactions  amongst  the  lattice  points,  and  in  an 
external  field. Further  simplification  is  achieved  by  assuming  that  there 
are  only  pair  interactions,  and  furthermore,  we  may  restrict  ourselves 
to  only  nearest  neighbor  interactions . In  this  case  a  simple  adsorption 
isotherm  can  be  obtained  from  the  Bragg-Williams  mean  field 
approximation  [$] 

o  * 

Xp(0)[g(<O] 

9  -  _  (13) 

a  •}  d/z 

i  +  Ap(°) tg(o) ] 

where  p’(O)  is  the  wall  contact  probability  of  the  SWP,  and 

g(<x)  is  the  contact  value  of  the  pair  distribution  function  for  the 

fluid, and  v  is  the  number  of  next  nearest  neighbors  of  the  lattice. 

* 

In  terms  of  the  Ising  field  B, magnetic  moment  _pt,and  interaction 
parameter  w  we  have 


Consider  now  the  case  in  which  a  hard  sphere  fluid  is  in  contact  with 


the  [111]  face  of  an  closed  packed  (hexagonal  or  face  centered  cubic) 
crystal  of  spherical  atoms  of  exatly  the  sa~a  diameter  as  those  of  the 
fluid. This  is  the  case  in  which  the  adsorbate  is  commensurate  with  the 
adsorbing  surface .Then, we  have  explicit  expressions  for  both  these 
quantities  [9,10] 


2. 


j<?(0)- jo(  i+iy/a-y) 

(16) 

g(o)-  (l+t?/2)i/(l-rj)C< 

(17) 

7  ‘ 

Replacing  these  expressions  into  (13)  we  get 


i  .  v  S 

A^(l+29)/(l-^)[(l+y2)/(l-rj)  ] 

9  “ _  (18) 

1  +  >p(l+2n)/(l-ip  [  (l+y2)/(l-^  )* 


Figure  1  shows  the  dependence  of  the  fraction  of  occupied  sites  9  with 
the  bulk  density  p  for  various  values  of  the  stickiness  )^.For  low  values 
of  ^  it  is  clear  that  there  is  a  high  degree  of  cooperativity  in  the 
adsorption, since ,  at  a  given  bulk  density  p  , the  fraction  of  occupied 
sites  9  changes  suddenly,  just  as  the  magnetization  of  the  Ising  model 
for  a  given  value  of  the  applied  magnetic  field. For  high  values  of  the 
parameter  ^  this  discontinuity  disappears,  and  the  isotherm  has  more 
the  apect  of  the  original  Langmuir  isotherm.  In  fact, for  low  densities 


we  get 
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